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LINEAR STABILITY OF COMPRESSIBLE VORTEX SHEETS IN 
TWO-DIMENSIONAL ELASTODYNAMICS 

ROBIN MING CHEN, JILONG HU, AND DEHUA WANG 


Abstract. The linear stability of rectilinear compressible vortex sheets is studied for 
two-dimensional isentropic elastic flows. This problem has a free boundary and the 
boundary is characteristic. A necessary and sufficient condition is obtained for the linear 
stability of the rectilinear vortex sheets. More precisely, it is shown that, besides the 
stable supersonic zone, the elasticity exerts an additional stable subsonic zone. A new 
feature for elastic flow is found that there is a class of states in the interior of subsonic 
zone where the stability of such states is weaker than the stability of other states in the 
sense that there is an extra loss of tangential derivatives with respect to the source terms. 
This is a new feature which Euler flows do not possess. One of the difficulties for the 
elastic flow is that the non-differentiable points of the eigenvalues may coincide with the 
roots of the Lopatinskii determinant. As a result, the Kreiss symmetrization cannot be 
applied directly. Instead, we perform an upper triangularization of the system to separate 
only the outgoing modes at all points in the frequency space, so that an exact estimate 
of the outgoing modes can be obtained. Moreover, all the outgoing modes are shown 
to be zero due to the L^-regularity of solutions. The estimates for the incoming modes 
can be derived directly from the Lopatinskii determinant. This new approach avoids the 
lengthy computation and estimates for the outgoing modes when Kreiss symmetrization 
is applied. This method can also be applied to the Euler flows and MHD flows. 
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1. Introduction 


Vortex sheets are interfaces between two incompressible or compressible flows passing 
along each other. They arise in a broad range of physical problems in fluid mechanics, 
aerodynamics, oceanography and astrophysical plasma. Some typical examples include the 
sharp interface between two parallel shear flows, and vortex flows where the vortices are 
concentrated within a thin layer. In particular, for compressible flows, vortex sheets are 
fundamental waves which play an important role in the study of general entropy solutions 
to multi-dimensional hyperbolic systems of conservation laws. Analyzing the existence 
and stability of compressible vortex sheets may shed light on the understanding of fluid 
dynamics and the behavior of entropy solutions. 

In this paper, we are concerned with the vortex sheets in the following two-dimensional 
compressible inviscid flow in elastodynamics ( [TSl |22l [3T] ): 


pt + div(pu) = 0, 

(/Ou)i -|- div(/9U (g) u) -|- Vp = div(/?FF^), 
Ft -|- u • VF = VuF, 


( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 


where p stands for the density, u = (u,n) G the velocity, F = {Fij) G the 

deformation gradient, and p the pressure where p = p{p) is a smooth strictly increasing 
function in (0, -|-oo). This model arises from the viscoelastic fluids with negligible viscosity, 
and describes the features of both fluids and solids. 

Before going into further discussions of the elastodynamics model, we would like to first 
review some known results in the study of the compressible vortex sheets. For vortex 
sheets in the compressible Enler flow, the study of stability in the linear regime dates 
back to 1950’s by Miles [43( I44j and Fejer-Miles |19] . It is a classical result that, in two 
or three spatial dimensions with Mach number M < y/2, the vortex sheets are violently 
unstable (c.f. [m |49]). These instabilities are the analogue of the Kelvin-Helmholtz 
instability for incompressible fluids. For the theory of the incompressible vortex sheets, 
we refer the readers to (UEKSgisollMlIiaiinilMllSH] and the references therein. In a 
series of papers, Artola-Majda ElElIll investigated the interaction of the vortex sheets 
and highly oscillatory waves in the two-dimensional compressible Euler flow, indicating 
the global-in-time nonlinear instability of vortex sheets for Mach number M > y/2. These 
instability results show that one can not expect the global existence of the vortex sheets 
in multi-dimensional spaces. In the pioneer works [Ml IT], Coulombel and Secchi used 
a micro-local analysis and Nash-Moser method to establish the linear and local-in-time 
nonlinear stability in two dimensions with March number M > y/2. In their setup, the 
initial data Is a small perturbation of a rectilinear (i.e. piecewise constant) vortex sheet 
and their notion of linear stability is in a similar sense to that of shock waves by Majda 
[39( 140] and Coulombel [MKH]. Their results imply the local-in-time existence of such 
vortex sheets. For two-dimensional non-isentropic Euler flows, Morando-Trebeschi [45| 
obtained a weaker linear stability of the vortex sheets, that is, there is some additional 
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loss of derivatives in their estimates. Moreover, there have been a lot of studies on the 
vortex sheets in the multi-dimensional steady flows; see p m Ea EZ] and the references 
therein. 

In complex fluids, the situation becomes more complicated. Ruan-Wang-Weng-Zhu [38] 
considered the linear stability of compressible vortex sheets in two-dimensional inviscid 
liquid-gas two-phase flows. They showed that the linear stability, compared with the 
Euler flow, may be weaker when the Mach number satisfies a certain condition. For the 
two-dimensional isentropic magnetohydrodynamics (MHD), Wang-Yu |56] showed that 
the magnetic fields will lower the critical Mach number and exert a small subsonic zone 
where some linear stability holds for the rectilinear current-vortex sheets. Moreover, for 
the three-dimensional compressible MHD, Blokhin-Trakhinin [7], Trakhinin |5IIE11E3| 
and Chen-Wang mm adopted a different symmetrization approach to obtain certain 
linear and nonlinear stability, and the existence of the current-vortex sheets. The results of 
MHD indicate the stabilization effects of the magnetic helds on the current-vortex sheets. 

For the viscoelastic fluids, there have been extensive studies on various aspects from the 
modeling and analysis point of view HSlIMllSllETlIlelllT], as well as on their applications 
HHEHIEHI. In the two important examples, namely, the shear flows and the vortex flows, 
both numerical experiments and theoretical analysis indicate that the viscoelasticity plays 
a stabilization role (see laiziiEsiiisi and the references therein). Moreover, vortex sheets 
in viscoelastic fluids have also been discussed by Huilgol [271 [28| , where the author consid¬ 
ered the Rayleigh problem in viscoelastic fluids and showed by constructing an example 
that the unsteady shearing motions can lead to vortex sheets in some viscoelastic liquids. 
On the other hand, Hu-Wang |26] showed the formation of singularity and the breakdown 
of classical solutions to system (iri])-(lL3|) for certain initial data. These results motivate 
us to consider the inviscid elastic fluids and investigate the stabilization effects of the elas¬ 
ticity on the vortex sheets. Specihcally, we obtain the necessary and sufficient conditions 
for the linear stability of rectilinear vortex sheets in a two-dimensional compressible isen¬ 
tropic inviscid elastic fluid in the sense of m by analyzing the Lopatinskii determinant 
of the linearized boundary value problem. 

As in the aforementioned works of Euler flow, two-phase flow and MHD flow, a common 
challenge of the vortex sheet problem is that the system has a free boundary, and the free 
boundary is characteristic. The characteristic boundary leads to some loss of control on 
the trace of the characteristic parts of the solutions HSIESIIII]. Moreover, the Kreiss- 
Lopatinskii condition does not hold uniformly, which implies some loss of the tangential 
derivatives in the estimates of the solutions in terms of the sources on the right side of the 
linearized problem naiiaEsi. 

In addition to the above difficulties, for the elastic flow, the appearance of the elasticity 
leads to a more complicated distribution of the roots in the Lopatinskii determinant. More 
precisely, the non-differentiable points of the eigenvalues may coincide with the roots of the 
Lopatinskii determinant, which does not happen for Euler flow (c.f. HSj). In the standard 
arguments da ESI ED, one needs to construct the Kreiss symmetrizer and work along with 
the Lopatinskii determinant at each point in the frequency space. At the non-differentiable 
points of the eigenvalues, the usual Kreiss symmetrizer requires no loss of derivatives with 
respect to the source term; but the degeneracy at the roots of the Lopatinskii determinant 
implies the loss of the tangential derivatives, which leads to some serious obstacles to apply 
the Kreiss symmetrizer argument. To overcome this difficulty, we deal with this problem 
in a different way. More precisely, instead of using Kreiss symmetrization to separate every 
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single mode, we first perform an upper triangularization of the system to separate only the 
outgoing modes from the system at all points in the frequency space. Then we establish 
an exact estimate on the outgoing modes from the equation, and use the L^-regularity 
of solutions to conclude that the outgoing modes in the homogeneous system are zero. 
As a result, we only need to estimate the incoming modes, which can be derived directly 
from the Lopantiskii determinant. Therefore combining the estimates for the outgoing and 
incoming modes we achieve the linear stability. In conclusion, we find that with the added 
elasticity, a new stability region can be generated in the subsonic zone, indicating the 
stabilization effect of elasticity as expected. We also find that within the stable subsonic 
region there exist a class of states where the stability of such states is weaker than the 
stability of other states in the sense that there is an extra loss of tangential derivatives, due 
to the fact that the Lopatinskii determinant exhibits higher order of degeneracy at such 
states. This is a new feature which Euler flows do not possess. We further remark that 
our upper triangularization simply tells us that the outgoing modes are all zero. Thus we 
can avoid the lengthy computation and estimates for the outgoing modes when the Kreiss 
symmetrization is applied, and hence the arguments are greatly simplihed. This upper 
triangularization method can also be applied to the Euler and MHD flows, as illustrated 
briefly in Section [71 

The result in this paper concerns with the linear stability with constant coefficients 
of the elastic vortex sheets. The natural next step is to consider the linear stability with 
variable coefficients, and more interestingly, the nonlinear stability of elastic vortex sheets, 
which would be a crucial tool for proving the local existence. These issues will be addressed 
in a forthcoming paper. 

The rest of the paper is organized as follows. In Section O we formulate the nonlinear 
problem of the vortex sheets, fix the boundary, linearize the system around a given trivial 
solution, state our main theorem and sketch the idea of the proof. In Section O we perform 
some preliminary reduction and eliminate the front, which leads us to a system of ordinary 
differential equations. In Section (H we consider the normal mode analysis and separate 
the modes. In Section (U we investigate the Lopatinskii determinant and introduce the 
estimates near the roots of the Lopatinksii determinant. In Section [6l we combine the 
separation of the modes and the estimates from the Lopatinskii determinant to obtain the 
energy estimate, and hence finish the proof of our main theorem. In the last section, we 
apply our method to the vortex sheets in Euler and MHD flows. 


2. Eormulation of the Problems and Main Results 

In this section, we derive the governing dynamics of the vortex sheets from the elas- 
todynamical equations (ll.ll) - (ll.3p . linearize them around a rectilinear vortex sheet, and 
state our main theorem. 


2.1. Governing equations for vortex sheets. To obtain the equations of the vortex 
sheet, we first reformulate the system (dUl-dOl) into a divergence form. Notice that the 
intrinsic property div(/9F^) = 0 holds at any time throughout the flow if it is satisfied 
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initially [25]. We can rewrite the system as the following form: 

pt + div(pu) = 0, 

< (/3u)t + div(/9U (g) u) + Vp — div(/9FF''") = 0, (2.1) 

_ {p¥j)t + div(/9Fj (g) u - u (g) pFj) =0, 

where Fj is the jth column of the deformation gradient F = [Fij), and i,j = 1,2. In 
column-wise components, the intrinsic property actually means 


div(/3Fj) = 0 for j = 1,2. (2.2) 

Consider U{t, xi,X 2 ) = (p, u, F)(t, xi, X 2 ) to be a solution to the system (12.11) which 
is smooth on each side of a smooth hypersurface T = {x2 = ijj{t,xi)}. We denote by 
12 = {—di'ijj, 1) the unit normal vector on F and 


U{t, Xi,X2) 


xi,X 2 ), when X 2 >'ip{t,xi), 

U~{t, xi,X 2 ), when X 2 <'ip{t,xi), 


where = (/9^,u^,F^). It follows that U satishes the Rankine-Hugoniot conditions at 
each point of F: 


dMp] - [pu ■iy]=0, (2.3a) 

dt-iplpu] - [(/9U • z/)u] - \p\v -F [pFF^u] = 0, (2.3b) 

dMp^j] - [(u • ’^)pFj] + [{pFj ■ z^)u] = 0, (2.3c) 

[pF, ■u] = 0, (2.3d) 


where [/] denotes the jump of the function / across the hypersuface F. Here we have used 
(|2.2p to derive the last jump condition. Now, we denote rrii, = p{ijjt — u - v). By (I2.3ap . we 
have [rriu] = 0. Thus combining (I2.3dj) and (I2.3cp . we can obtain 


'mu[Fj] + {pFj ■ i^)[u] = 0. 


(2.4) 


Since we consider the contact discontinuity, we assume dttp = u+ • v = u“ • v. By ()2.3ap 
and ()2.4p . we obtain [p] = 0, rriy = 0, [Fj • z/] = 0 and Fj-v = 0, for j = 1, 2. Therefore for 
a vortex sheet in the elastic flow, the jump conditions ()2.3ap - ()2.3dp become 

p'^ = p~ , V’t = u”*" • ^ = and F'^ ■ v = FJ ■ u = 0, j = 1,2 on F. (2.5) 

We now introduce the transformations ^^(t, xi, X 2 ) to straighten the free boundary F as 
follows. We first consider the class of functions <I>(f, xi, X 2 ) such that d 2 ^(t, xi, X 2 ) > k > 0 
and <h(t,xi,0) = ijj{t,xi). Then we define 

U# = (p^,u^,F^)(t,xi,X 2 ) = (p,u,F)(f,xi,4>(t,xi,±X2)) 

for X 2 > 0. For simplicity of notation, we drop the ^ index in the rest of the paper. 
Dehne ‘h^(t,xi,X 2 ) = ^(t, xi, ±X 2 ). To obtain a simpler formulation of the problem, we 
can choose <1> satisfying 

= 0 , 


for any X 2 > 0, which is inspired by [161120] . With this change of variables, we can rewrite 
the equations (|l.ip - (ll.3l) in the following form: 

dtU^ + Ai{U^)diU^ + [A2{U^) - dt<h^I - di<h^A3iU^)] d2U^ = 0 ( 2 . 6 ) 
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for X 2 > 0 with fixed boundary X 2 = 0, where 
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0 

0 

0 \ 


p 

V 

0 

-Fii 

0 

—Fi2 

0 


0 

0 
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0 
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0 
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0 
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0 

0 

-Fii 

0 
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0 
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0 
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0 

0 

0 

V 
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\0 

0 

— Fi2 

0 

0 

0 

V / 



i “ 

0 

p 

0 

0 

0 

0 \ 


0 

u 

0 

—F 21 

0 

—F 22 

0 


F 

p 

0 

u 

0 

—F 21 

0 

—F 22 

A2{U) = 

0 

—F 21 

0 

u 

0 

0 

0 


0 

0 

—F 21 

0 

u 

0 

0 


0 

—F 22 

0 

0 

0 

u 

0 


Vo 

0 

—F 22 

0 

0 

0 

u / 


Moreover, from (12.5p . we obtain the boundary condition 

' (u+ — v~)di'ip — (tt+ — u~) = 0 
= 0 

< F+^dii) - ^21 = 0 , at X 2 = 0, (2.7) 

(F+ - - {F+ - F2-2) = 0 

F+SiV’ - F+ = 0 

^p+-p- = Q 


where ^ at X 2 = 0. 

It is easy to see that the above system (I2.6h - (j2.7j) admits rectilinear solutions. Moreover, 
all the rectilinear solutions can be transformed under the Galilean transformation and 
change of the scale of measurement to the following form: 


( 





F 

0 


0 

pr 

, U-:= 


0 


0 

F[2 


FI 2 



\ 0 y 


^^{t,Xi,X2) = ±X2, 


where the constants p, , v\ F[]^, FI 2 and ^(2 satisfy 


( 2 . 8 ) 


F + v^ = F[i + F[^ = FI 2 + F (2 = 0 and F > 0, Fl^ / 0. 


2.2. Linearized problem. Now we linearize the system (|2.6p - (l2.7l) aronnd the above 
constant states (j2.8p . Denote = (/9^,u^,F^) = and to 











STABILITY OF VORTEX SHEETS IN ELASTODYNAMICS 


7 


be the small perturbations of the constant solution, and consider the following linearized 
problem: 

dtU^ + Ai{U^)d\U^ ± A2{U^)d2U^ = 0, 

in X 2 > 0, with the boundary condition at X 2 = 0: 


' (u*' — v'')diip — (t(+ — ii ) = 0, 
dt(p + v^diif — = 0, 

{F{,-Fl,)diip-iF+-F,-,) = 0, 
< F[^di^-F+=0, 
{F{,-Fl,)di^-{F+-F,-,) = 0, 
— -F 22 “ 

- P~) = 0 , 


where = (p at X 2 = 0. In short, we have 


L'U = 0, if X2 > 0, 
B{U,ip) = 0, if X2 = 0, 


where 




B{fi,p) = 


U-J ' V 0 Ai{U-) 

/ {v''— v‘)dip — {u'^ — u~) \ 
dtp + v^dip — 

{FA-Fl,)d,p-{F+-FF,) 

F[,dip - F+ 

iFl2-Fl2)d,p-iF+-FF2) 

F[,dip - F+ 

V {P^ - P ) J 


u+ 

(j- 


+ 


'A2{U+) 

0 


0 

-A2{U- 


u+ 

(j- 


(2.9) 


Next we symmetrize the system (I2.9h . Consider the following change of variables, 

T 0\ (U+ 

0 T [u- 


W = 


( 2 . 10 ) 


where T is a matrix of the following form: 

/ 0 


T = 


2p 

0 

0 

0 

\ 0 


0 

J_ 

2c 

J_ 

2c 

0 


0 0 0 0 \ 
0 0 0 0 
0 0 0 0 
10 0 0 


0 0 10 0 
0 0 0 1 0 
0 0001 / 


and c = ^Jp'{p) is the sound speed of the constant solution. For convenience, we denote 
the components of the new variable by 

IF = (IFi,IF2,IF3,-,1^14)^, 
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and the tangential, normal, characteristic and non-characteristic parts of W by 
lytan ^ We, 1^11, W"l3)^, 

= (1^2,1^3, We, W 7 , Wg, hbio, 1^12, Wu}^, 

= (1^1, W 4 , We, We, Wj, TVs, Wii, f^i2, hhis, Wuf, 

W^<^ = iW2,W3,Wg,Wio)^. 

With the above change of variables, we multiply the system (12.91) by a symmetrizer 
A = diag{l, 2 c 2 , 2 c 2 , 1 , 1 , 1 , 1 , 1 , 2 c 2 , 2 c 2 , 1 , 1 , 1 , 1 }. 

Thus we can obtain the following equation 


where 



'AqL'U 

= £W 

= AodtW + AidiW + A2d2W = 0, 



B{U,^) 

= B(W-,<p) = 

MW'^ + 6 ( 1 

\di^J 

= 0 

? 


-4.1 = 1 

"A[ 0 ^ 

A 

with 








/ 

-C2 

c2 

Tpr,l 

-^11 

0 

-^12 

0 

\ 


-C2 


0 

0 

cnii 

0 

-rF"'’^ 


c2 

0 2c^v'^’^ 

0 

-cF^’^ 

cr^ll 

0 

771 

“'^-^12 

A]:^ = 

_ 7?^’^ 

-^11 

0 

0 

vr,l 

0 

0 

0 


0 

- 

C2*11 

C2*ii 

0 

vr,l 

0 

0 



0 

0 

0 

0 

vr,l 

0 


\ 0 

-cF^’^ - 

-cF^’^ 

0 

0 

0 



A = diag{0, —2c^, 2c^, 0,0,0,0,0,2c^, —2c^, 0,0,0,0}, 




c —c 

0 0 

c c 

0 o\ 


/o 

2v^ \ 



—c —c 

0 0 

0 0 

0 0 


1 




-1 1 

0 0 

1 - 

1 0 0 


0 

0 


M = 

0 0 

-1 0 

0 0 

1 0 

, b = 

0 

to 



0 0 

-1 0 

0 0 

0 0 


0 

F[i 



0 0 

0 -1 

0 0 

0 1 


0 

2F[2 



0 0 

0 -1 

0 0 

0 0 ) 


VO 

FI 2 ) 



( 2 . 12 ) 


We now introduce some weighted Sobolev spaces needed for the main theorem. Define 
:= {u{t,xi) G D'(M2) : e-^^u{t,xi) G 
:= {vit,xi,X2) G : e-^^v{t,xi,X2) G 

for s G M, 7 > 1, equipped with the norms 


respectively, where 


:= ||e 


\Hi 


n := l|e '^^vWhs 




= {(t,a:i,a:2) G : X2 > 0}. 


( 2 . 13 ) 









STABILITY OF VORTEX SHEETS IN ELASTODYNAMICS 


9 


In our stability analysis, the norm that is frequently used is 


I^IIa7 = ( 2 , 


JR 2 


where u{^) is the Fourier transform of u. Note that by letting u := e we have that 
II^IIh®(r2 ) and Hfills.^ are equivalent, denoted by ||rt| 




— m 


s,7‘ 


Now we can define the space equipped with the norm 


7 

H-oo 


llL2(ffi) 


f+co 

/ \\v{-,X2) 

Jo 




dX2. 


Again we have 




r+oo 

' := / 

Jo 


\v{-,X2)\\l^dX2. 


Note that || • ||o ,7 is actually the usual norm on L^(]R^) and 

L2(M3 ). 


l0,7 


is the usual norm on 


2.3. Main result and methodology. With the above notations and spaces, our stability 
result can be stated as follows. 

Theorem 2.1. (1) If the particular solution defined by (12. Sp satisfies 

{v’-f > 2c" + (Ffi)" + (Ffj)", or 

, r,2 , ,r,r 7 , , r-r ,2 . , , r.2 , {(Hi)'" + (H 2 )") (2c" + (Hi)'" + (H 2 )") (^'l^) 

(v ) < (F„) + (Fij) . but (v ) -4((Fri)" + (Fy2 + c")-’ 

then there is a positive constant C such that for all 7 > 1, VF € and ip € 

the following estimate holds: 


(2.15) 


L2(H0) + ll^“U2=0|li2(R2) + ||<F||^1 

(2) If the particular solution defined by (12.Sp satisfies 


^ 4((F7)2 + (Ff 2 ) 2 +c 2 ) 


or 


(2.16) 


then there is a positive constant C such that for all 7 > 1, VF € and p € 

the following estimate holds: 


h{HO) + ll^“U2=0|li2(]R2) + II‘F||^1(R2) 

< c (T|||£if|||J.,„., + 

(3) If the particular solution defined by p2.8h satisfies 

{Fr = {Ff,f + {Ff,f + 2cfi 


(2.17) 


(2.18) 
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then there is a positive constant C such that for all 7 > 1, VK G and (/? G 

the following estimate holds: 

7|||W^IIIl 2(H0) + ll^’^U2=o|li2(R2) + ||¥^||^1(R2) 

(4) If the particular solution defined by (j2.8p satisfies 

< 2c2 + (Ff,? + {Ff,)\ 

the constant vortex sheets (j2.8p is linearly unstable, in the sense that the Lopatinskii con¬ 
dition is violated. 

Remark 2.1. The above theorem gives all the situations (l)-(3) where the linear stability 
holds. On the other hand, for the remaining situation (4), we can show that the linearized 
problem is unstable due to the failure of the Lopatinskii condition. The details will be 
discussed in Section [5l Therefore this theorem actually gives us a sufficient and necessary 
condition on the stability of the linearized problem. 

Remark 2.2. (i) From the second condition of (I2.14p we observe that the elasticity term F 
produces a stable subsonic region, and hence provides a stabilization effect on the vortex 
sheets. It is easily seen that in absence of the elasticity F = 0, our results recover the 
stability theory for the vortex sheets in the Euler flow. 

(ii) Another feature of the elastic flow which is different from the case in the Euler flow, 
can be inferred from the first condition in (I2.16p . It follows from (I2.16P that there is a 
class of states in the interior of the subsonic region where the stability holds in a weaker 
sense, i.e. there is a two-order loss of derivatives in the estimates (12.171) . 

(iii) The second condition in (12.161) and the condition (j2.18l) give the transition states 
across which the stability property changes. This property is similar to the Euler and 
MHD flows (c.f. [151 [56]). At these two classes of states the stability holds in a weaker 
sense. 

Now we sketch the idea of the proof of the main theorem. We follow the standard 
argument (c.f. m) to first remove the source term from the equations, eliminate the 
wave front ip from the resulting system, and then single out the non-characteristic part 
of the unknown W to arrive at a reduced system in the Eourier space of the form: 

< ( 2 . 21 ) 
JW^%,=o = h, 

where A is a 4 x 4 block diagonal matrix and /3 is a 2 x 4 matrix (the explicit forms are 
given in (14.2p and (IS.lip i. It turns out that all the desired estimates in Theorem 12.11 can 
be achieved by estimating W^‘^\x 2 =o- 

By counting the number of boundary conditions one can only hope to control at most 
two components of W'^^\x 2 =o- Thus in order to obtain the full estimate on 13,2=0 one 
has to utilize the differential equation in (I2.2ip . The conventional way to do so is to use the 
Kreiss symmetrization. Roughly speaking, this argument is to first find a 4 x 4 Hermitian 


(2.19) 


( 2 . 20 ) 
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matrix r and a number C > 0 such that r + 0(3*(3 and rA are both positive definite. The 
positivity of r + C(3* (3 implies 


< > 1^2=0+ C|/ip > K 




X2=0 


( 2 . 22 ) 


for some k > 0, where < •, • > is the Euclidean inner product in C^. Therefore if one can 
further choose r such that < > \x 2 =o < 0i then VT“| 3 ; 2 =o can be controlled by 

h. To check < > 1^2=0 < 0, we multiply the differential equations in (I2.2ip by 

r*lT“ [r* is the Hermitian transpose of r) and integrate with respect to X 2 from 0 to oo 
to obtain 

1 

~ 2 

Thus a sufficient condition for stability is that rA is positive definite. However at the non- 
differetiable point of the eigenvalues of A, the usual symmetrizer r can not make r 3-0/3*13 
positive definite for any C > 0 when the Lopatinskii determinant is zero there, and hence 
the Kreiss symmetrization method seems hard to apply. It turns out that this situation can 
happen for the vortex sheets in elastodynamics. More precisely, if the particular solution 
defined by (12.8p satisfies 

Kf = t{Wi)" + (Frj)" + c") (2.24) 

with c < v'3((Ffi)2 + (T{ 2 ) 2 ), some non-differentiable points of A will coincide with some 
roots of the Lopatinskii determinant (c.f. Remark 15.II) . and hence the Kreiss symmetriza¬ 
tion cannot be applied directly. We will develop some new idea to overcome this difficulty 
which we now describe as follows. 

Our approach is to first perform an upper triangularization of A to obtain a closed 
differential system of the two components of which are not in the stable subspace of 
A (we refer to these two components as outgoing modes and the other two components as 
incoming modes). This way the differential system (I2.21jl is transformed to 


W--^rW^‘^>\x,=o = 


f 


< W^^,rAW^'' > dx2. 


(2.23) 


d I 

d^ 


(G * \ / W*" \ 
\^0 H j \wout j ’ 


(2.25) 


where W*” and are both two-dimensional vectors corresponding to the incoming and 
outgoing modes of IT“, G and H are both 2x2 matrices. The closed differential system 
of is 

dX2 

This upper triangularization ensures that all the eigenvalues of H have positive real parts 
(in fact for our case, H can actually be a diagonal matrix whose diagonal entries both 
have positive real parts). Hence an exact estimates of can be obtained; furthermore 
by the L^-regularity of IT“, one has that = 0 for X 2 > 0. Therefore to obtain the 

full estimates for IT“|a; 2 = 0 ) h remains to estimate the other two components of IT“|a; 2 = 0 ) 
i.e. lT*”|a; 2 = 0 ) which are in the stable subspace of A. Instead of integrating the first two 
rows of (j2.25p to derive the estimates of W'‘'^\x 2 =o-, we only use the boundary conditions 
in (I2.2ip . In fact, from W°'^^\x 2 =o = 0, the boundary conditions are reduced to 

= h, at X 2 = 0, 


(2.26) 
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where P is a 2 x 2 matrix whose determinant is exactly the Lopatinskii determinant. Thus, 
if the Lopatinskii determinant is not zero, P is invertible and the matrix norm of P~^ can 
be estimated. Therefore fT*”'|x 2 =o = p-^h, and VL*”U 2 =o is controlled by h. This together 
with the fact that 3 , 2=0 = 0 lead to the estimates of W'^‘^\x 2 =o and hence complete 

the proof of stability. We further point out that our new approach can be applied to 
other fluid models including the Euler and MHD flows, which will be illustrated in the 
last section. 


Before we conclude this section, we introduce the following change of unknowns in order 
to simplify the notations in our proof of Theorem 12.11 

W = ex.p{—'yt)W, if = exp(— 

with 7 > 1. Denote two new operators C"^ and by 

= e-^^CW = -iAoW + AodtW + Aid^W + W, 

= MW^ + b ■ 

We note that and ||e“'>'*u||s ,7 are equivalent to the norms and ||u||j?s 

respectively. Then we have the following equivalent formulation of Theorem 12.11 In the 
later context, we aim to prove the following theorem. 


Theorem 2 . 2 . (1) If the particular solution defined by (12. 8 p satisfies (|2.14p . there is a 
positive constant C such that for all 7 > 1, VE G and G the following 

estimate holds 

2 


7 


__ 

2 

__ 

IT 

+ 

VT"U,=o 


0,7 


0,7 


+ 


2 

1,7 




C^W 


H —2 

1,7 7 


^ ^^(^"13,2=0,^?)"" 


1:7 


(2.27) 


(2) If the particular solution defined by (12. 8 p satisfies (I2.16p Then there is a positive 
constant C such that for all 7 > 1, W G and G the following estimate 

holds 


7 


W 


0,7 

< c 


+ 




0 : 2=0 


0,7 


+ 


2 

1>7 


/ 1 

ow 

2 1 


H- 7 

V7^ 


to 

4 ; 


B^(W“l3:2=0,^) 


2,7 


(2.28) 


(3) If the particular solution defined by (|2.8p satisfies (I2.18P Then there is a positive 
constant C such that for all 7 > 1, IT G i7^(Mi]_) and G 77^(]R))_), the following estimate 
holds 


7 


IT 


0,7 
< C 


+ 


IT" 


1 0 : 2=0 


0,7 


2 

1,7 


C^W 


+ 


1 


3,7 T 


B^{W^\x2=o,t) 


3,7 


(2.29) 


(4) If the particular solution defined by (12.81) satisfies 

< {Pf < 2 c 2 + + (P[2)^ 


(2.30) 




































































STABILITY OF VORTEX SHEETS IN ELASTODYNAMICS 


13 


the constant vortex sheets (j2.8l) is linearly unstable, in the sense that the Lopatinskii con¬ 
dition is violated. 


3 . Decomposition of the System and Elimination of the Front 

In this section, as a first step for the proof of Theorem 12.21 fand hence Theorem 12.ip we 
perform some preliminary redactions and eliminate the front ip from the linearized system. 
Consider the following inhomogeneous differential system for W and p on M5]_: 

r c^w_= i if X2 > 0, 

\B'y{W^,p)=g, ifx2 = 0, 

where / and g are given sources. 


3.1. Decomposition of the system. Due to the linearity of the system mu , we can 
decompose it into two subsystems as follows. First we consider the following problem for 
V with the homogeneous boundary conditions: 


where 


CW = f, if X2 > 0, 
MiD° = 0, ifx2 = 0, 


Ml 


/o 1 0 0 
0 0 10 
0 0 0 1 
0 0 0 0 
0 0 0 0 
\0 0 0 0 


0 0 0 o\ 
0 0 0 0 
0 0 0 0 
10 0 0 
0 0 10 
0 0 0 1 / 


(3.2) 


From the classical hyperbolic theory [B], the boundary condition is strictly dissipative, and 
thus (j3.2l) admits a solution satisfying following estimates: 


2 C, 

Ss-i 


lir"k=ollI„ < ^ 


2 

A;,7’ 


for any integer A: > 0. 

Then we consider the second problem for the difference W = W — V. In fact W 
satisfies the following homogeneous differential equations with inhomogeneous boundary 
conditions: 


aw = 0, if X2 > 0, 

13^{W^,p)=g:=g-M_W', if X 2 = 0 . 


Remark 3.1. Here we are slightly abusing notation, since W was previously defined to be 
the transformed perturbation of the rectilinear vortex sheets, c.f. (I2.11h . From now on we 
will consider IT as a solution to (|3.3I) . 


Multiplying the equations in system (13.3p by W and integrating, we have 
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Thus it suffices to derive the following estimate on W: 


C 




|2 

lfe ,7 


(3.4) 


to obtain all the estimates in Theorem 12.21 where k will be determined accordingly. 

Now we perform the Fourier transform to system (13.3p with respect to tangential vari¬ 
ables {t,xi) and denote the variables in the frequency space by (<5, r/). Let t = ^ + i6. 
Then W, the Fourier transform of W, satisfies the following system: 


{tAq + ir]Ai)W + ^ 2 ^— = 0 , if X 2 > 0 , 
dX2 


(3.5) 


, 6 (r, 77 )^ + MW'^ = g, 


if X2 = 0 , 


where 

/ \ 

T + iv’^rj 
0 

6(r, g) = h- = 2iF[^g 

iF[,g 

V 

To utilize the homogeneity structure of the system, we define the hemisphere 
S = {(r, r/) : \t\^ + = 1 and > 0 } 


in the whole frequency space 11 := {(r, 77 ) : r € C, 77 £ M, |rp -(- 77 ^/ 0 , > 0}. It 

is easily seen that 11 = {s • (r, 77 ) : s > 0, (r, 77 ) G S}. We will carry out our argument 
on the hemisphere S and use homogeneity property to extend it to the whole frequency 
space n. 


3.2. Elimination of the front. An important observation in the vortex sheet system is 
that the wave front (/? is only involved in the boundary conditions. Thus we can estimate 
the wave front by using the ellipticity of the boundary conditions of (|3.5I) in the sense of 
the following lemma. 


Lemma 3.1. There is a C°° map Q : 11 —^ GLi{C) which is homogeneous of degree 0, 
such that 


Q{T,g)b{T, 77 ) 


/ 0 \ 

0 

6 »(r, 77 ) 

0 

0 

0 

V 0 / 


V(t,77) Gn, 


where 9 is homogeneous of degree 1 and 


min \9{t, g)\ > 0. 

(u'?)GS 
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Proof. We only sketch the proof of the lemma. The idea is to consider the map 


/ 0 

0 

1 

0 

0 

0 

0 \ 

r + iv^T] 

—2iv^r] 

0 

0 

0 

0 

0 

—2iv^r] 

f — iv’^rj 

0 

— 


-2iFl^g 


-Ffi 

0 

0 

v'^ 

0 

0 

0 

-Ffi 

0 

0 

0 

2v^ 

0 

0 

-FI 2 

0 

0 

0 

0 


0 

\ -F[2 

0 

0 

0 

0 

0 

2v^ / 


on S and extend it by homogeneity of degree 0 to the whole frequency space. Then the 
lemma can be proved by a direct computation. □ 


Now we multiply the boundary conditions in (13.5p from the left by Q and obtain the 
new boundary conditions: 

Qbif + QMW'^ = Qfj, at X 2 = 0. (3.6) 


With this choice of Q, we have QM = 



1 

0 

0 

1 

-1 

0 

0 \ 

—c(t — iv’^g) 

—c(t — iv^rf) 

0 

0 

c(r + iv’^r/) 

c(t + iv^rf) 

0 

0 

—c(f — 3iv'~ri) 

—c{t — 3iv'~r]) 

SiF[,g 

3 iF[^V 

—2civ^r] 

—2civ^r] 

- 2 iFf,g 

- 2 iF[^g 

cF[, 

cF^i 

-v^ 

0 

-cF^i 

-cFf, 


0 

cFf, 


-2v^ 

0 


-cFf, 

0 

0 

cFh 

cFf 2 

0 

—v'' 

-CFI2 

-cF[2 

0 

v'' 

CFI2 

cFf 2 

0 

-2v'^ 

-CFI2 

-CFI2 

0 

0 / 


on S which is homogeneous of degree 0. Denote the third row of the above matrix by i. 
Hence the third equation of the new boundary conditions (13.6p is 

eij, r])(p + i{e, r])W^\x^=o = b*{e, r])g, on S. 


Notice that i and b* are homogeneous of degree 0. From the compactness of S and 
continuity of ^ and 6*, we know that they are bounded on H. By Lemma l3.II and a direct 
integration of the above equation, we have 

||y^||?,,<C(||W"U,=o||§ + Nlg). (3.7) 

Moreover, the last four equations in the new boundary conditions (13.61) are 



cFfi 

—v^ 

0 

-cF{, 

-cF[i 

v’" 

0] 

cFfi 

cFfi 

-2v^ 

0 

-cFf, 

-cF[i 

0 

0 

CFI 2 

^Fl2 

0 

—v'" 

-CFI 2 

-cF{2 

0 


VF{^ 

CFI 2 

0 

-2v^ 

-CFI 2 

-cF{2 

0 

0/ 


for all (r, rj) G H, where Q is the matrix consisting of the last four rows of Q. Isolating 
the third, forth, seventh and eighth columns of the left matrix in the above equations, we 
can obtain 


/ 

—u'’ 

0 

u’’ 

0\ 


/IF 5 \ 



-2v^ 

0 

0 

0 


W7 



0 


0 

u'’ 


Wu 


V 

0 

-2v'^ 

0 

V 


[WuJ 



/Ff, 

F[i 

-F[i 

-F[A 

F{i 

Ffi 

-F[i 

-F[i 

F{2 

Fr2 

-F{2 

-Fr2 

[f{2 

Fr2 

-F{2 

-FI2I 




'a: 2 = 0 ) 
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on n. Since now the left matrix above is constant and invertible, we can obtain 

2 


1x2=0 < C 




X2=0 


+ l?l' 


(3.8) 


Recall the definition of in (I2.10p . We know from ()3.8p that an estimate of || 0 : 2 = 0 IIo 
leads to the estimate of ||kR“U 2 =o|lo> hence the estimate of ||(/9||f from (j3.7p . Therefore 

by (I3.4p we know that it is sufficient to obtain the estimate of ||hk“‘^U 2 =o|lo- 

Now consider the first two rows in the new boundary conditions ()3.6p at X 2 = 0 and the 
equations of (j3.5p for X 2 > 0. We arrived at the following system where the wave front (p 
is eliminated: 


Qj VV 

{tAq + ii]Ai)W + A 2 -;— = 0, 
dX2 


(3.9) 

(3.10) 


where h is the product of the first two rows of Q and and 

\—c{t — iv’^rf) —c{t — iv’^rf) c(r + iv'^rf) c(r + iv'^rf) 

on S and is extended to all of 11 by homogeneity of degree 0. Note that from the homo¬ 
geneity of Q we have the following bounds 


(3.11) 


\h\^<C\g\\ 


(3.12) 


for some positive constant C which is independent of The rest of this paper is 

devoted to the estimate of ||iy’^'^ 13 , 2 = 0 ||o from the system (I3.9I) - (|3.10I) . 


4. Normal Mode Analysis 

In this section, we consider the normal modes and want to separate the outgoing modes 
from the system microlocally. This separation gives us the exact estimates on the outgoing 
modes. Moreover we will show in Section [6] that these outgoing modes are all zero. 

4.1. Normal modes. To obtain an estimate of ||IT“U 2 =o|Ioj to derive a closed 

differential system of IT“. For this purpose, we single out the following ten algebraic 
equations from (13.9p : 

(r + iv^r])Wi - ic^r]W2 + - iF^2VWe = 0, 

- iFiirjWi + (r + iv'"ri)W4, = 0, 

- icF{^r]W2 - icFl^gW^ + (r + iv''g)W5 = 0, 

- iFl2f]Wi + (r + iv'^vi)WQ = 0, 

- icF{2VW2 - icFl2gW^ + (r + iv^j])W7 = 0, 

(r + iv^r])Ws — ic^rjWg + ic^rjWiQ — iF^gWu — fT(2??VFi3 = 0, 

- iFli^Ws + (r + iv''r])Wn = 0, 

- icF^^rjWg - icFl^rjWio -|- (r -|- iv''r])Wi 2 = 0, 

- iFl2'nWs + (r + iv’-T])Wi3 = 0, 
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— “ icF'uV^io + (t + iv^r])Wi4^ = 0 . 


Recall definition (|2.11l) . we know that the above equations can be explicitly solved in terms 
of = {W 2 , 1 ^ 3 , as the following: 

ic^ri{T + iv^r]) 


IRi = 
W 4 = 
1^6 = 
tR 8 = 

IRii = 
1 Ri3 = 


(t 

+ iv^pfi 

+ {{FhY + 

(F[2: 

j2)^2V 



—YFf^^r? 



(t 

+ iv^pY 

+ {{FfiY + 

(FI 2 : 

j2)^2l 



-c^Fl2ri^ 



(t 

+ iv^pY 

+ {{FfiY + 

WiV. 

)2)^2l 


k 

Yr]{T + iv^rj) 


f 

{'T 

+ iv^pY 

+ {{FliY + 

(FI 2 ) 

|2) r/2 ^ 



—YF[.^r? 




^ + iv^rjf 

■ + ((Fhr + 

(F 




-YFl^r/^ 




{W 2 - W 3 ), 


icFl^r] 

T + iv^T] 

icFl^rj 

T + iv^rj 


{W 2 + W 3 ), 


(W 2 + W 3 ), 


{Wg - IRio), 


(Wg - IRio), Wu = (^9 + ^ 10 ), 

r + iv'-rj 


(Wg - IRio), fRl4 = 


icFl^r] 


(Wg + lRio). 


(r + + (^^2)^) ^ T + iv'^T]' 

Then using the differential equations in (|3.9I) , we obtain the following differential equations 
for only; 

d 


= AW^ 


where 


dx2 
/ n'’ —m 


A = 


0 

0 


o\ 

0 


0 0 —n m 

y 0 0 —m^ v} j 


(4.1) 


(4.2) 


and 


r,l 

n ’ = 


2 (r + + ((Ffi')^ + (T[^^)^) 


r I ^ 

m ’ = - 


2 c(r + iv^’^rj) 
(r + iv'^'^r])rf‘ 


+ 


(t + iv^’^ri)rj 


’ (r + zu^>^r/)2 + (^(F[{^)2 + tf 

({Fit? + (Fit?) 


(r + + ((T [/)2 + (Fit?) r ?2 2 c(r + iv^h^) 


From the classical hyperbolic theory (see, for example, 0), it follows that one of the 
key issues in the estimates of ||kF°'^|x 2 =o||o is to bound the components of on the 
stable subspace of A. For this reason, we first derive the following lemma of Hersh-type 
|23j on the explicit description of the stable subspace of ^ on S. The proof can be done 
by a direct computation, and hence we omit it. 


Lemma 4.1. For {T,rj) £ S and 91r > 0, ^ defined in (14.2p admits four eigenvalues ± 0 ;'’ 
and Fio^, where and are negative. Moreover, the following dispersion relations 




















18 


R. M. CHEN, J. HU, AND D. WANG 


hold: 


{u'-f = {n’-f - (m'-f = 4 [(T + + ((Hi)" + (U 2 )") »"] + ri\ 


c 

= {n^f - (mf = ^ [(r + iv^V? + {{F^f + 


+ rf. 


(4.3) 


The eigenvectors of oj'^, uj^ and —oj^ take the following forms respectively: 

EL = (a", 6 ", 0,0)^, E\ = (a^c^0,0)^, 

E^_ = ( 0 , 0 , 6 ^a')^, E^^ = ( 0 , 0 ,c^ay, 

where 


(4.4) 




If,I = _ a;"’V’', 

.r,l\2 


a 


T,l 


= {t + iF’^p) (r + iv^’^vf + ((Fltr + (Fl^r ) h 


:,r,l\2\ 2 


Both of and oj^ can he extended continuously to all points (r, r/) G S such that 9^r = 0, 
so can EL and E^- The two vectors EL and EL remain linearly independent for all 
(r,? 7 ) G S. 

By the definition of A, (14.31) actually holds on 11. Moreover, from (14.41) . we can see 
that A can not be smoothly diagonalized near the points (r, 77 ) G S satisfying one of the 

following: = 0 , or = 0 , or r = ±iv^g, 01 t = i ± ??, 

because EL and E^ (or EL and E^^) become parallel at those points. In fact, (r, ry) G S 

are poles of A when r = Eiv^p or t = i ± ^ {Ff^)'^ + ( 1 ^ 12 )^^ Therefore instead 

of looking for a diagonalization of A, we perform an upper triangularization of A, which 
we refer to as separation of modes. 

4.2. Separation of modes. Since the eigenbasis of A may degenerate at some points on 
S, we need to treat A microlocally. This means that for each point (tojT/q) G S, we will 
separate the outgoing modes of A from the system (14.1|) in some open neighborhood V C S 
of that point (to,77o)- Here we refer to the outgoing modes of A as all the components 
of which do not belong to the stable subspace of A. By using this separation, we 
can show later in Section El that for every (To,r/o) £ T, these outgoing modes in fact 
vanish in V H {(r, 77 ) : 3f?r > 0}. Thus from the compactness of S, we can extract a 
finite covering of S to show that the outgoing modes of A vanish in the entire 

S n {(r, 77) : > 0}. 

The following proposition guarantees that the separation of modes can always be achieved 
for all points on S. 

Proposition 4.1. For defined in Lemma\4.1\ we have 


for all (r, 77) G S. 


(r + iv^’^r])u^’^ - c -v"^) 7^0 


Proof. The key idea of the proof is to examine the signs of the real and imaginary parts 
of For this purpose, we consider the general situation where {x + iy)"^ = p + iq for x, 
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y, p, q € M and x < 0. By a direct computation, we can express x and y in terms of p and 
q as 


X = —1 


lp+ + ? 


y = -sgn(g) 



(4.5) 


for all (p, q) € M^\{p < 0 , g = 0 }. 

Next we apply the above relations to at the point (r, 77 ) G B with Sftr > 0. Let 
^r,l _ ^r,l + where x^’*, y^’^, p'’’^, q^’^ G M. From the definition 

(j4.3ll of a;'’’* we know that x^'’* < 0 and 

f = Ml + Lh). (4.7) 

With this setup it is easy to see from (I4.5I) that when {p^'\q^’^) 0 {p < 0,q = 0} and 
6 + 7 ^ 0 the sign of p'"’* is opposite to the sign of 5 + v^’^rj respectively. Here we recall 

the fact 7 = 3ftr > 0. On the other hand, when (p^’^ ^ {p < 0,q = 0}, (I4.5p fails to 

hold. More precisely, those points correspond exactly to the case when 

7 = 0 , 5 + / 0 and p'’’^ < 0 . 

Therefore they all lie on the boundary of B. By Lemma I4.ll the values of at the 
boundary of B are defined as the continuous limits of the interior values of This 
way, we can still determine the signs of x^’^ and by continuity. Thus by the continuous 
extension of along the path where the ratio of S and y is fixed, the sign of is opposite 
to the sign of (5 + v^'^y respectively at those exceptional points (p’’’^, g'’’*) G {p < 0, g = 0}. 
Hence we conclude that 


if 5 + v^'^y / 0 , the sign of p^’^ 


is opposite to the sign of (5 + v'^'^y respectively. 


(4.8) 


Now we return to the proof of the proposition. We will only prove that (r + iv'^y)uj'^ — 
c ((w^)^ — p^) / 0. The case that c ((rc^)^ — — {T + iv^y)uj^ / 0 can be treated similarly. 

Assume on the contrary that 


(r + iv^y)uj^ — c — y^) = 0 . 


(4.9) 


If r + iv^y = 0, the above equation becomes (cu'’)^ — = 0. Combining this with ()4.3p . 

we obtain ((Tf^)^ + (F^lj)^) y^ = 0, and hence p = 0, which in turn implies r = 0. This 
contradicts the fact that (r, p) G B. 

Thus we can assume r + iv^y 7 ^ 0. From this it follows that 


c ((o;'')^ — p^) 
(r + iv^y) 


1 

c 


(r + iv'^y) + 


+ {Fl^f) y^ 

T + iv'^y 


(4.10) 


If 7 (i.e. 5Rr) is positive, it is easy to check that the real part of the right hand side of the 
above formula is positive, which contradicts the definition that < 0 . 

Thus we only need to check the situation when 7 = 0 . In this case, we have r + iv^y = 
i((5 + v^y) ^ 0. By ()4.10p . we know that = 0, and therefore = 0 and p^ < 0. We 
further claim that p^ 7 ^ 0. Otherwise if p^ = 0, from (14.3h and the fact that g*" = 0 we 
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have = 0. Then from (j4.9l) we must have r] = 0. By (I4.6p it follows that 5 = 0, and 
hence r = r/ = 0 , which contradicts the fact that (r, ?]) G S. 

Therefore we only need to consider (r, r]) £ T, when r + iv'^rj 7 ^ 0, 7 = 0 and < 0. 
This immediately implies that 6 + v'^rj 7 ^ 0. From (14.81) we know that the sign of 
is opposite to that of 5 + v'^rj. However by (14.1011 and the fact that = 0, 

c {6 + v'"r]) 

Since p"^ < 0, from (14.Op we know that {6 + v^r])'^ — > 0- Thus the sign 

of is the same as the sign of 5 + v^rj, which is a contradiction. Hence we conclude 
that (r + iv^r])uj^ — c ((w'’)^ — 7 ^) never vanishes in S, which completes the proof of the 
proposition. □ 


With this proposition, we can show that E'p!' do not vanish at any point on S. Because 
if = 0, we have = 0 and {rf'^ — = 0. By a direct computation we 

have that 7 ^ 0. Then = 0 implies that rrf’^ = 0. From the definition of rrF’^ it 

follows that 


c (r + iv'^’^ri)rf‘ 

^ (r + iv^'^ 7 ])^ + + (72)^) 


((7i7 + (F'Jr) 7 

2 c(r + iv’'’^r]) 


Together with = 0 and (14.3p . we have 

(r + - c - 772 ^ = 0 , 

which contradicts Proposition 14.11 

The non-degeneracy of allows us to construct the following transformation matrix 


T = {EL,F^,Ei,F^} 


and similarly. 


= 


F' = 


in a neighborhood of (to,t/o) S S with 

(0,1,0,0)T, 
(1,0,0,0)T, 

(0,0,1,0)T, 
(0,0,0,1)T, 


if 74 0 at (to,t/o)) 
if (n^ - 74 0 at (ro,7o), 


if m}a^ 74 0 at (to,?7o), 
if (n* — Lij^)a^ 74 0 at (to,77o)- 


Obviously, from the above argument, for any point (to,Po) S S, there is an open neigh¬ 
borhood V of (to,Po) where T is continuously invertible. Then we can obtain 


T-^AT = 


on V where A is given in (j4.2p and 

1 


2’'-' = 


a' 


m 


r,l 




0 0 \ 
0 -w'’ 0 0 

0 0 Lo'' z’' 

000 -uy 


if 74 0 at {to,po), 

if (re^’* — 74 0 at {to,pq). 


(4.11) 
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Remark 4.1. By the definition of rrf'\ and we know that may not be well 
defined at the poles of A which are all located on the boundary of B. However, the 
estimates of W'^‘^\x 2 =o only involves the interior points of B. Hence it sufficies to obtain a 
uniform esitamtes of W’^‘^\x 2 =o hi the interior of B, which corresponds to B H {3ftr > 0}. 


5. Lopatinskii Determinant 

In this section, we want to estimate the components of IH"‘^| 3 ; 2 =o in the stable subspace 
of A through the boundary conditions, which requires us to investigate the invertibility of 
the matrix (3{E^, EL)- This can be done by computing det{(3{E^, E^_)), which is known 
as the Lopatinskii determinant. By a direct computation, we can simplify the Lopatinskii 
determinant to be 

A = det{/3{E^, E’_)) = c^{t + iv'^r]){T + iv^rj) ((r + — c ((w'’)^ — rj'^)) 

(e: 

(c - 772) - (r + + a;'), 

from which we see that the Lopatinskii determinant A can vanish at multiple places in 
B, which indicates that one can not expect the uniform Lopatinskii condition to hold. 
In the following lemma we detail the root distribution of A, which is important for later 
discussion. 


Lemma 5.1. The roots of the Lopatinskii determinant A are distributed in the following 
ways. 

Case 1. If v'^ > + (^{ 2 ) 2 , then all the roots are simple and on the 

boundary ofT,. The Lopatinskii condition holds. More precisely, the roots are {T,r]) G B 
such that 

1 . T = diiv^'q, 


2 . r = 0, 

3. T = dziVir], 

where + (Ef^)^ + (^[ 2 )^ + c2 - + 4+ (^[ 2)2 + c2) (v^)T 


Case 2. If 0 < + (-^ 12 )^ but ^ ^ 

then all roots are also simple and on the boundary of B. 
holds. More precisely, the roots are (r, ry) G B such that 


{iFIir+{F[2)^){2A + iFIi)^ + iFl2)^) 
The Lopatinskii condition still 


1 . T = diiv^'q, 

2 . T = EiVig. 


r. n r, r / ( 2 c 2 + (Ei-B2+(Ef2)2) i. a 

Case 3. It v = \ - - /A then all roots are on the boundary 

V 4((F{'B2 + (Ff2)2+c2) 

ofT,. The Lopatinskii condition still holds. More precisely, the roots are {T,rj) G B such 
that 

1. T = dziv’^rj (double roots). 


Case 4. //u'’ = y^2c2 + (T{^)2 + (F|2)2, then all roots are on the boundary o/B. The 
Lopatinskii condition still holds. More precisely, the roots are (r, ry) G B such that 

1 . T = Eiv'^rj (simple roots), 

2 . T = 0 (triple root). 
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Case 5. If , then all roots are on the boundary of S. The 

Lopatinskii condition still holds. More precisely, the roots are (r, r/) € S such that 

1. T = ±iv^r] (simple roots), 

2. T = 0 (double root). 


Case 6. If ^ + {F 12 Y < + (^[ 2 )'^, then some roots are in 

the interior of S, and hence the Lopatinskii condition fails. 

Moreover for the roots above, if the degree of a root (r, rj) ^ T,, with r = idr] for some 
real number d, is k, then we have A = (r — idr])^h{T,r]) for some continuous h{T,r]) 
satisfying h{T, r?) / 0 near the root. 

Remark 5.1. It is easily seen that if 

= j{Wl)" + (U2)" + <!") (5.2) 

with c < y^3 + (.^^ 12 )^) (this corresponds to the second case of the Lemma ISTTI) . then 

one of the non-differentiable points of A, namely t = i (.^n)^ + {F 12 Y + 

coincides with the root r = iv'^r] of the Lopatinskii determinant. As we have pointed out 
in the Introduction, this is a new phenomenon which does not appear in the compressible 
Euler flow. 


Proof. The proof of the above lemma depends on a detailed analysis on each factor of the 
Lopatinskii determinant. Firstly, the factors (r + iv^'^r])uj^'^ — c are exactly 

the expression in Proposition 14.11 Thus we know they are never zero. 

Secondly, we consider the factors r + iv^’^rj. Obviously r = —iv'^’^rj are the only simple 
roots to r + iv'^’^rj = 0 , respectively. 

Thirdly, we assume that 

U;'^J - 7^2 = 0 . (5.3) 


If r/ = 0, we have w'’ = oj^ = — ^ ,which means 7 ^ 0, and hence uj^lo^ — rj"^ 7 ^ 0. 
Moreover, if d + v’^'^rj = 0, for example d + v^rj = 0, then from (j4.3p uj^ is real and negative. 
However, since 77 7 ^ 0, we know that d + 7 ^ 0, which implies 7 ^ 0. Thus can 

not be a real number, which violates (15.3p . Therefore (j5.3p can not happen for r/ = 0 or 
d + = 0. 

This leads us to focus only on the points where g 0 and d + v^’^'g 7 ^ 0. Introduce the 
following two variables. 


From ()5.3I) we have 


V = 


= 


CO 


r,l 


ig ig 

— 1 , and hence (H’') 2 (H ^)2 = 


1. By the (14.3p we know 


(5.4) 


(Si’')2 = I[(r + o f - (F[,f - (F[,fj - 1, 
(fi')" = Il(V + ff - (Fl,f - (Fl,fl - L 

Hence we have 


[(F +77^)2 -(F[,)2 


- c2][(P + v^ - - (Fl2f - C^] = 


(5.5) 

(5.6) 
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which leads to the following equation for V^: 

_ 2 ^ ^ ^ ^ 2 ) ^2 ^ ^ r 4 _ ^ ((^ r ^)2 ^ ^ ^ 2 ) (^.)2 

+ ((Ffi)' - (F[2ff + 2c2 ((F[i)2 - = 0. 

Using the quadratic formula, the two roots of the above equation are 


= (v^f + (F[,f + (Fl^f + - ^c4 + 4((Ffi)2 + (F[2)2 + c2)(u’')2, (5.7) 

U| = + (F[,f + (F[2)2 + c2 + ^c4 + 4((F[i)2 + (F[2)2+c2)(u>')2. (5.8) 


We claim that the points (r, ??) G S with r = ±iV 2 r] are not the roots of (15.31) . Without 
loss of generality, we can assume V2 is positive. By a direct computation, we have V 2 +v^’'' > 
\J(? F (7^11)^ + (7^12)^ and -U2 + 1^'’’^ < + (F[^)2 + (^{2)^. Now, by (143)) . (14771) and 

(BSD, have opposite signs to d + respectively. If r = iV2??j we have 7 = 3f?r = 0 

and 5 + v'^’^rj = (V2 + u'’’*)^. Therefore w'’ and iv* are purely imaginary, and 


from which we deduce that 


= 




G M, 


sgn 


= - 


sgn 


6 + 

77 


= -sgn ( U2 + = -1. 


Therefore 7^ ~1 and (|5.3p is not satisfied. Similarly, we can show that (r, 7) G S 

with T = —iV 2 T] are also not the roots of (15.3p . 

Now we focus on Vi- Obviously, by (15.7p . we have 


U2 > 0, if > ^2c2 + (F[i)2 + (i7r^)2 or 0 < < ./{FJ^ViFT^-, 

V,^ = 0, if u*- = ^2 c 2 + (Ffi)2 + (F[2)2 or = ^(F[,)2 + (Ff^)^; 

U2 < 0, if y/ 2 c 2 + (F[i)2 + (F[2)2 > > ^(F[i )2 + (F[2)2. 


If y^2c2 + ( 37^)2 + ( 3 ( 2 )^ > v'^ > {F[^)‘^ + ( 372 )^, by (15.41) . we obtain that r = izUi?] 

are real. Thus (5 = 0, but 7 / 0 and Sir / 0. This implies that 6 + v^’''t] ^ 0. By (14.6p 
and p4.7p . we know that and = —q^ / 0. Using (14.5p . we can have 

and y*' = —y\ i.e. is the complex conjugate of uK Then uj'^uj^ > 0, which implies that 
r = FiViT] are the roots of (15.3p . This way we are able to find a root (r, 77 ) with > 0, 
which violates the Lopatinskii condition, and hence the vortex sheets are unstable. This 
proves Case 6 in the lemma. 

For the rest cases we will consider > 0 and when taking the square root we always 
use the positive branch, that is Vi > 0 . 

If v'^ > + ( 37^)2 + (^{ 2 )^, we have r = FiVip is purely imaginary. Without loss 

of generality, we only consider the case when r = iVip. Then = 0, but <5 7 ^ 0 and 
7 / 7 ^ 0. By a direct computation, we can obtain 


Ui+U^’' >^(3fi)2 + (F[2)2+c2. 


( 5 . 9 ) 
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Thus 6 + v^’^r] = {Vi + v^’^)r] / 0. By (|4.3p and (j5.9p . we have that are both real 

and negative. This means are purely imaginary, and from (14.81) the signs of are 
opposite to those of (5 + respectively. Hence 

sgn = —sgn ^(5 + v'"r]){S + v^r])j = —sgn ^(Vi + v'^){Vi + = 1 . 

Therefore oj^oj^ > 0, and (r, ?/) G S with r = iViij are roots of (15.3p . The other case when 
r = —iViT] can be treated exactly the same way. 

If v'" < + (T’l’ 2 )^, we also have that r = dziVii] is purely imaginary. Similarly 

as before we only treat the case r = iVir], Then 3f?r = 0, but <5 7 ^ 0 and ?7 / 0. Now, we 
have 


Vi+v 


rl 


< 


0Ffi)2 + (F[2)2 + c2. 


(5.10) 


By (14.3p and (I5.10p . we have are both real and positive. Thus are both real 

and negative, which implies > 0. Hence (r, r/) G S with r = TiViry are roots of (15.3p . 

Then we want to show that under the condition n'’ > y/2c^ + (T^i)2 + (^{’ 2)2 or 0 < 
n'’ < sj (T|'^)2 + (^{ 2 )^ the roots to (15.31) are simple. Since (15.31) does not admit a root at 
r/ = 0, the points (r, 77 ) G S satisfying = 0 are not the roots of — 77 ^ = 0 From 
(j4.3l) . a;^’* are analytic near the points where do not vanish. We can differentiate (15.5p 
and (j5.6p with respect to H at F = Hi to obtain 


Thus 


dV y=y^ H^’T2 


+ 1 ) 

dF 


y=Ui 


c 2 H^O' 


Plugging in (15.51) and (15.61) . we obtain 


d(0^0' + 1) 


dV 

Using P5.7I) . we have 


2Fi (Fi 2 - (F)2 - (Ffi)2 - {Fl^r - c2) 


y=yi 




d{n^n^ + 1) 


dF 


/O. 


y=yi 


Hence we have proved (r, rj) £ T, with r = ±iVir] are all simple roots of (15.3p provided > 
y^2c2 + (Tfi)2 + (^{’ 2 )^ or 0 < . More precisely, near r = ±iFi 77 , 

we have — 7/2 = (r ± iVirj)h^{T,r]) for some continuous /i^(r, 77 ) / 0 respectively. 

If = y^2c2 + (T{'^)2 + (^{ 2 )^, we obtain r = ± 7 Fi 77 = 0. In this case = d = 0 
but 77 / 0 . By p4.6l) . we have = —(? < 0. Together with the fact that = 0, 
d + 7 ^ 0 and (14.8p . we infer that are purely imaginary and have 

opposite signs to respectively. This implies > 0. Hence (r, 77 ) G S with r = 0 
are roots of p5.3p . 

If = y^(F{j^)2 + (^{ 2 )^, we also obtain r = ± 7 Fi 77 = 0. Then = d = 0 but 77 7 ^ 0. 
By p4.6l) and p4.7p . we have > 0 and = 0. This implies are both real and 

negative. Thus ui^ui^ > 0. Hence (r, 77 ) G S with r = 0 are roots of (15.3p . 
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Now we want to check the multiplicity of the root when r = 0 under the condition 
n'’ = y^2c^ + + (-^ 12 )^ or . Similarly as in the previous case, 

we obtain the following first derivative of + 1: 


+ 1) _ 2y (y2 _ _ (^rj2 _ ^^r^)2 _ ^2) 

dV ~ 

Further differentiation yields the following second derivative: 

+ 1) _ 6F2 _ ((^'■)2 + ^ ^2^ ^2V {V^ - K)" - - {Fl^f - ^)y 


[dVY 

Thus 




+ 1 ) 


dV 


d?{FFFl^ + 1 ) 


= 0 , 


V=Vi 


{dVf 


/O. 


y=yi 


Hence (r, t]) ^ T, with r = 0 are all double roots of (15.3p if = ^J2c^ + + (^{ 2 )^ 

or n'’ = + (^{ 2 )^. We can also conclude that = T^h{T,r]) for some 

continuous h{T, ??) / 0 near r = 0. 

Finally we are left to consider the last factor of ()5.1[) 

uj^ + J = 0. (5.11) 


It is obvious from the (14.3p that if > 0 then < 0 and < 0, and hence 
Lo^ + Lo^ ^ 0. Thus we focus on the case when 3ftr = 0. By (j4.7p . we have = 0. From 
(j5.1ip and using the definition (j5.4p . we have which implies = pK Using 

we have 

2v^6p = 2v^5p. 


This implies Sr] = 0. If rj = 0, we obtain 5 = 1. Thus p"^ = p^ = —^ < 0 and S + = 

1 > 0. From (|4.8p . we obtain < 0, which contradicts (|5.11l) . Hence it must be 

(5 = 0. Since {t,p) € S, we have p = ±:^. In this case, 




If + (^{ 2 )^ + — {v^y > 0 , we obtain that are both real and negative, which 

contradicts (j5.1ip . Otherwise if + (^{ 2 )^ + cP — (v^y = 0, we have lo^’^ = 0, which 

means (0,±4^) are roots of (15.111) . On the other hand, this situation belongs to Case 6 
which we have already concluded the emergence of the instability. Therefore in the sequel, 
we only consider the case + (^( 2 )^ + P — (v^y < 0. In this case are purely 

imaginary. Since 3f?r = (5 = 0 and p = ±:^, we have S + v'^p = —{S + v^p) p 0. Again a 
use of ()4.8[) implies that the signs of are opposite to those of S + v^'^p, which implies 
Sjcj’' = —Hence w'’ + w* = 0, and (0,±^) € S are roots of (j5.1ip . 

Next we turn to the multiplicity of these roots. From (|4.3p we know uj^ and w* can not 
vanish simultaneously. Thus + oj^ is analytic near these roots. From ()5.4I) . ()5.5I) and 
dM]), we have 

dnr,i Y F 


dV 
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This implies 

d{n^ + n^) _v + v^ v + v^ 

Since in this case = —oj^ / 0 , we have fi’' = 7 ^ 0 at ( 0 , ±;j^). At these roots, we 

have 

dV ^ 

Hence, (0, ±^) are all simple roots to (j5.11jl . if + (-^{ 2 )^ + 0 ^ — < 0. Therefore 

we have = Th{T, rf) for some continuous /i(t, ?/) / 0 near r = 0 . 

To summarize, we have derived all possible roots (r, rf) of the Lopatinskii determinant, 
namely, 

r = —r = izHir/ or r = 0, (5-12) 

where we have assumed that > 0. In general, some of the roots may coincide, and we 
have already discussed the possibility that Hi = 0 in the study of (15.31) . Now we are left 
to check whether v'^ = Hi when Hi > 0. By a direct computation, = Hi if and only if 


V = 


' (( F [,)2 + ( 2 c 2 + ( F [,)2 + ( F [ 2 )^) 

4((F[,)2 + (Ff2)2 + c2) 


This way, we can identify the root distribution (j5.12j) with the six cases in the lemma: 
Hi^ > 0 and Hi / u'’ Case 1, Case 2, 

Hi^ > 0 and Hi = Case 3, 


Hf = 0 Case 4, Case 5, 
Hi^ < 0 =)> Case 6 . 


Putting all the above together we finish the proof of the lemma. □ 

With this lemma, we can obtain the following estimates on the stable subspace of A 
near the roots of the Lopatinskii determinant. 


Lemma 5.2. Let (To,r/o) £ be a root of the Lopatinskii determinant A. Then there is 
a neighborhood V o/(to,? 7 o) which does not contain any other roots of A and a constant 
kq, such that for V(r, 77 ) G V and VZ“ G 
If > v'2c2 + (F[i)2 + (Ff2)2, 


If0<v^< V(Ffi )2 + (Ff 2 ) 2 , but 7 ^ 

\P{E^,E^)Z-\‘^ > ko7^|Z-|2 


{(FLr + {F[2)^){2c^ + iFhr + {F[2)^) 
4(Hfi)"+(iT2)"+c") ^ 


Ifv^ 


I {{FL)^+iFLr){2c^ + {FL)^ + {F[,r) ^ /,^rV2 , 717r N 

V A{{F[,r+{F[,r+c^) +1^12) 


\I3{E^,E^)Z-\^>ko7^\Z-\^ 
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If = V2c2 + 
when tq = hiv^rjo, 

> ao7^|Z-|2 

when Tq = 0 , 

Ifv^ = V(F[,)2 + {Fl^Y, 
when Tq = ±iv'^7]Q, 

\P{E^,E^)Z-\‘^ > ko7^|Z-|2 

when Tq = 0 , 


(3{E^,E^) = 


dll 

d2i 


di2 

d22 


Proof. We denote the elements in the Lopatinskii matrix by 

-a'’ + ¥ a} - 6' 

c(r — iv'^rf){al + 6'’) c(r + iv^r]){a^ + 6^)^ 

Since every element of (5{E’E^) is continuous, we know that if there is an element of 
I3{E^,E^) which is nonzero at (to, 7 o)) then there is an open neighborhood V of (ro, 7 o) 
such that I3{E’E^) can be transformed to a diagonal matrix in V, that is 


Pfi{EL,Ei)Q=(^^ 


0 

A 


(5.13) 


For some P, Q continuously invertible in V. For example, if dn 0, then we have the 
following identity 


( 1/dii 

y—d2i/dii 


fj P{E\E 



-di2\ 

dll ) 


1 0 
0 A 


Now we claim that there is always an element in fi{E'', E^) which is not zero for all (r, if) 
in S. First we consider 


dll 

di2 


[(r + iv^rjf + ((F[i )2 + (^[ 2 )^) 7 ^] [(r + iv^r])u}^ 
(r + iv'-r]f + + (^( 2 )^) rf\ [(r + iv'-rfjj 


c{{ui^? - V')] , 


From Proposition 14.11 (r + iv^ri)io^ — c{{io^)‘^ — if') and (r + iv^rj)uj^ — c((a;^)^ — 77^) are 
never zero. Thus dn only vanishes when r = —iv’^rj ± i\/{F[iy^I-{FJ^'il , and di2 only 
vanishes when r = ± + {Ff2)‘^f]- If du = di2 = 0, we can obtain r = 0 and 

= sj (Ff ^)2 + {E[2Y- Since (r, r/) £ S, we have 7 / 0 , which impiles that r ± iv^rj / 0 
and / 0 . From the expression of ^21 and ^22, we know that in this case ^21 / 0 and 
d22 / 0 . Hence for any (ro,7o)) there is an open neighborhood V of (ro,r7o) such that 
I 3 {E^, E^) can always be locally continuously transformed to diag{l. A}. Therefore by the 
continuity and boundness of dij, equation (j 5 . 13 l) implies that 

\P{EL,E’'_)Z-\‘^ > Amin(l, |A|2)|Z-|2, 


in V, where a > 0 depends on (ro, 7 o)- From the fact that there are only hnitely many 
roots of A we know that V can be chosen so that it only contains one root of A, which is 
(ro, 7 o)- Now combining the result in Lemma l5.11 we finish the proof. □ 
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Remark 5.2. For the points (to,? 7 o) where the Lopatinskii determinant is not zero, by the 
continuity of A, we can also obtain an open neighborhood V of (to,? 7 o) in which A 7 ^ 0 
such that 

\I3{EL,eL)Z-\^ > Ko\Z-\^, 

for all (t, 77 ) G V and Z~ G where kq is a positive constant depending on (To,r/o). 


6. Energy Estimates 


Now we want to combine the argument in the previous two sections and obtain the 
energy estimates. For a point (to,? 7 o) £ 5], shrinking the neighborhood if necessary, we 
obtain a new neighborhood V of (to,7?o) where we have the separation of modes of A (c.f. 
(14.111) 1 and the estimate of the Lopatinskii determinant (c.f. Lemma [5.2l and Remark l5.2p . 
We call such a point (to,?]o) a “generating point” of V. Notice that if A(to,? 7 o) / 0 then 
A 7 ^ 0 at every point of V. Repeating this process for all points on S forms a covering 
of S. Then by the compactness of the S, there is a finite subcovering of S with 

the corresponding generating points denoted by {(tj, 77 *)}^^. Obviously this subcovering 
contains all the neighborhoods V of (tq, 770 ) such that A(ro, 770 ) = 0. Then we can construct 

a partition of unity of S according to this finite subcovering, i.e., we can find Xi £ 

N 

for i = 1, • • • ,N such that ^ y? = 1 on S. 

Now we derive an energy estimate in each conic zone Ilj = {(r, 77 ) : s ■ {T,rj) G 
Vi, for some s > 0}. In each neighborhood Vi of {Ti,rji), we denote Ti the transfor¬ 
mation matrix of the separation of modes in this neighborhood and extend Xi Ti by 
homogeneity of degree 0 to the conic zone Ilj. Then we consider 

Z = XiT-^W'^'^ (6.1) 


= {t-^at,) Z. 


for all (r, 77 ) G Rj. The system that Z = (Zi, Z 2 , Z 3 , ^ 4 )"'' satisfies now becomes 

dZ 
dX2 

Recalling Remark l4.ll we know that it suffices to obtain uniform estimates of Z for IRr > 0. 
Therefore in the below we only consider IRr > 0. 

From (j4.11l) . the second and fourth equations are 


dZ2 

dX2 

dZ^ 

dX2 


-ufZ2 


— u} Zi^ 


for all (r, 77 ) G n* with > 0. By (14. 3 1) we have (r, 77 ) < 0 provided IRr > 0. 
Moreover, since VF(r, 77 , •) is in 1? and is continuous and bounded from above in Ilj, 
we know that Z{T,rj, •) is also in for every (r, 77 ) G Ilj. Hence, the above ODEs implies 
that 


Z 2 = 0 and Z 4 = 0, 


for all (r, 77 ) G Hj with IRr > 0. 


( 6 . 2 ) 
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Next for Z\ and Z 3 , from (16.11) and (16.211 we have 




XiW^‘^ = T,Z = {EL,E^_) 
for all (r, rf) S Ilj with IRr > 0. Then the boundary condition becomes 
X^h = x^PW^%,=o = P{EL,El) 


Zi 

Zz 


(6.3) 


3 : 2=0 


for all (r, ry) £ Ilj with IRr > 0. In the above equation, det (^(3{E^, EL)) is the Lopatinskii 
determinant A. From Remark 15.21 if det (/3(£^!L, F'L)) is not zero at we have 

\l3iEL,Ei)Z-\^>Ki\Z-\\ 

(r, 77 ) £ Vi and Z~ £ where k* is a positive constant depending on {Ti,rji). Since [3 is 
homogeneous of degree 0 , we have 

\l3{EL,EL)Z-\^>Ki\Z-\‘^, 

for all (r, r]) £ Ilj and Z~ £ . By (j6.3p , we have 

2 


(Z{\ 


[zJ 

X2=0 



(Zx\ 


CO 

to 

II 

0 


<—\hf 

for all (r, 77 ) £ Ilj with > 0 . 

If iTi,rji) is a simple root of A, by Lemma [521 we have 

\^{el,eL)z-\^ > Ka^\z-\\ 

for all (r, 77 ) £ Vi and Z~ £ M^. Since /3 is homogeneous of degree 0, we have 
(|r|2 + \P{EL,eL)Z-\‘^ > Ka^\Z-\\ 

for all (r, 77 ) £ Ilj and Z~ £ . By (j6.3p , we have 

KiY 

for all (r, 77 ) £ Ilj with IRr > 0. 

If {Ti,r]i) is a double root of A, from Lemma 15.21 we have 

\f3{EL,EL)Z-\^ > Ka^\Z-\\ 

for all (t, 77 ) £ Vj and Z~ £ M^. It follows from the homogeneity of /3 that 
(|t |2 + iv^W)^ \p{el,eL)z-\^ > Ka^\z-\\ 
for all (r, 77 ) £ Ilj and By (j6.3p . we have 

2 (M^ + (W )',„2 

< Xi ^- 4 -^ \h\ 

i^iT 

for all (r, 77 ) £ Ilj. 

Similarly if {Ti,rji) is a triple root of A, we have 

\P{EL,El)Z-\^>Ka^\Z-\\ 


(zA 


UJ 

3 : 2=0 


(6.4) 


(6.5) 


( 6 . 6 ) 
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for all (t, 77 ) G Vi and Z G Using homogeneity of /?, once again we have 

(| t |2 + > Ka^\Z-\\ 

for all (r, rf) G 11* and Z~ G . By (16.3p , we have 

_ ^ —sy— 

for all (t, 77 ) G 11*. 

Putting together (|6.2I) . (|6.4p - (l6.7l) we obtain the following estimate for Z in II*: 

2 ^ 2(l'^l^ + (^’')^^^)^i.|2 


(ZA 


[zs) 

X2=0 


(6.7) 


^L,=o ^X* 




2j 


\h\\ 


( 6 . 8 ) 


where j = 0 corresponds to the case when j5{E'^_,E^_) is invertible at (r*, 77 *), and j = 1,2, 3 
corresponds to the multiplicity of {Ti,T]i) as a root of A. 

Now we proceed to the proof of Theorem 12.21 

Proof of Theorem \2VA When (I2.14p holds, from Lemma l5.II we know either /3{E^, E^_) is 
invertible at (t*, 77 *) or (r*,r 7 *) is a simple root of A. Then from (j 6 . 8 p we can always have 

2 , 2(kl" + K)V),„2 




KiY 


\h\\ 


for all (r, 77 ) G 11* with iRr > 0 and all i = 1, • • • , A. From () 6 .ip . we have 


xl 




X2=0 


K-a 


Using the boundedness of E in II* and summing up the estimates over all the conic zones 
m we obtain 




3:2=0 


< C A!!±lpvi |ft|2 

T 


for all (r, 77 ) G II with > 0. Then integrating the above inequality with respect to (<5, 77 ) 
over and recalling (I3.12p we have 

2 C ^ 






which implies (I2.27p . 

The other two cases can be treated the same way. For the sake of completeness we 
provide the details. When (I2.16jl holds, Lemma [5T] indicates that (r*, 77 *) may be a nonzero 
point, simple root or double root of A. From (16.81) we always have 


I ^ 1 x 2=0 I 


KiT 


\h\\ 


for all (r, 77 ) G II* with iRr > 0 and all i = 1, • • • ,N. Thus from (16.ip . we have 

(|r|2 + (A)V)' 


xl 




3 : 2=0 


< Xi 






2 
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Patching up the estimates over all the conic zones {11./}^^ we obtain 

2 ^ 




\K 


for all (r, ?/) € 11 with > 0. Integration in (5, rj) over implies 




2 C*,, „2 

„S:^ll9lkr 


which proves (|2.28p . 

The last case when (j2.18p is satisfied, from Lemma fS.ll we know (rj, r]i) may be a nonzero 
point, simple root or triple root of A, and hence we have 

I VI |2 / .,2(1^1^ + 

I^L2=o| — Xi 6 1^1 ’ 

^11 


for all (r, rj) S 11 j with > 0 and all i = 1, • • • ,N. Converting into we have 

-—fi- \h\ . 


xl 




3:2=0 


Therefore 


VP” 


2 < (^( 1 ^ 1 ^ + 


Ki'y^ 

2 ^ 2\3 


<y\J 

for all (r, 77 ) G 11 with > 0, which shows that 


\hY 


VP” 




proving estimate (I2.29p . Thus we finish the proof of our main theorem. 


□ 


7. Applications to Other Models 

In this section, we provide two examples where our method of separation of modes can 
be applied. 


7.1. Euler. First, we consider the compressible vortex sheets in the following isentropic 
Euler flow (c.f. m) 

dtp + div(pu) = 0, 

dt{p\y) + div(/9u O u) + Vp = 0. 


We use the same change of variable $, and linearize the system around the rectilinear 
vortex sheets 







±3^2, 
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where v'' = —v^ > 0 and p are real numbers. We keep the same argument as in m up 
to Lemma 4.2. By m Lemma 4.2], we know that ca'" and have nonzero real parts, 
provided that > 0 on S, and the corresponding eigenvectors are 


EL{t, vj) = ^(r + - (r + 0,0^ , 

Eiir, p) = 1^0, 0, i(T + iv^pf + ^p^ - (r + iv^p)J, . 


By a direct computation, and note that ui^ and ui^ are the roots with negative real parts 
for > 0 on S, we have E'[_ ^ 0 and F}‘_ ^ 0. By the same argument as in Section [4.21 
for V(ro, po) G S, there is an open neighborhood V and a matrix T{t, p) which is invertible 
in V, such that 


TAT-^ 


/a;- 

r 

z 

0 

0 \ 

0 


0 

0 

0 

0 

J 


VO 

0 

0 

-ijj 


(7.1) 


Similarly as in our model, and may blow up to infinity at some points on the boundary 
of S. But since only the points (r, p) with > 0 will be involved in the energy estimates, 
the blow-up points again do not affect the stability. 

Next we turn our attention to the Lopatinskii determinant. By m Lemma 4.5], for any 
root (tq, ?7o) of the Lopatinskii determinant, there is an open neighborhood V of (tq, po) in 
S and a constant kq = ko{tq,po) > 0, such that following estimates hold for all (r, p) G V: 


MZ- G 


/ 3 ( T , r /) {eL{t,p),E^_{t,p)^ 



2 I ry — I 2 

> «o7 1^ I 


(7.2) 


For the points (to,t/o) which is not a root of the Lopatinskii determinant, we can still 
obtain an open neighborhood V of (ro,r7o) in S and a constant kq = ko{to,po) > 0, such 
that the following estimates hold for all (r, p) gV: 


VZ- G C^, 


PiT,p) {eL{t,p),E^_{t,p)^ Z 


> Ko|Z p 


(7.3) 


Thus by shrinking the open neighborhood V for each point {tq,pq) G T, if necessary, we 
can obtain the matrix T and the constant kq such that (|7.ip and (|7.2p or (j7..Sjl hold for all 
(r, p) G V. Hence we can apply the same argument as in Section [6] to recover the stability 
results as in m- 


7.2. MHD. The second example we consider is the following compressible vortex sheets 
in the following isentropic MHD (c.f. [56]): 

dtp + div(pu) = 0, 

dt{pu) + div(/9U (8 )u-H(8 )H)-hV(p-h ^|H|^) = 0, 

(9tH - V X (u X H) = 0. 
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We linearize the system around the rectilinear vortex sheets 


Ur 


! p\ 


f p\ 

0 


0 

^2 

,u^ = 

vi 

0 


0 

W 2 ) 


W 2 J 




±Xl, 


where the components of are constants satisfying 

v^^+vi = 0, v^^>0, \H^2\ = \H^ 


21 - 


Similar ly, keeping the same argument as in [56] up to Lemma 4.1, we obtain 

dX2 

/3W^%2=o = k 


where 


with 


A = 


(rf —rrf 0 0 ^ 

-rf 0 0 

0 0 —n^ im} 


/3(T,r?) = 


y 0 0 —m‘ n'' J 

-pX^ pX^ 


pX^ 


-P>? 


X{t + 1x2^) A(t + —X{T + iv2p) —A(r + ’ 


rrf' = — < — 




+ iv’/p) 


+ 


A I 2{T + iv2^p) 2 (^X‘^{t + iv'/pY + 


r,l ^ ) f , ■ r,l ^ , CA^ri^ 

n ’ = - pr + iv^ 77) + --+ 


2 (r + zug S) 2 ^A 2 (r + zz ;2 ^77)^ + j 


X = Jc^ + c\, CA = I P- 


By [561 Lemma 4.1], if (r, 77 ) E S with > 0, the eigenvalues uf and lJ have negative 
real parts and the corresponding eigenvectors are 

EL{t,p) = (a"77z",a"(n''-a;"),0,0)^, 

E^_{t,p) = (^0,0,a\n^ — uj’‘),a^m'‘'^ , 

where 

= (r + iv^’^'p) ^X^(t + iv^'^p)"^ + c^ca^p'^'^ ■ 

Then we want to examine whether E'^’^ can vanish on S. Setting E^ = 0, it is obvious 
that d^{rf — rrf — lo’’) = 0, which, by |56l Lemma 4.2], it is equivalent to the condition 
that a'' = 0 or (r + iv’^pY + c\p = 0. Note that rrf is undefined at a'’ = 0. Moreover 

lim rrf 7 ^ 0 . 
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On the other hand, a direct computation shows that when (r + iv^r])'^ + = 0, we have 

7 ^ 0. Thus 7 ^ 0 in S. Similarly we can show that E[_ 7 ^ 0 in S. Hence again we can 
apply the argument as in Section [4 . 21 to upper triangularize A in an open neighborhood V 
of each point (To,ryo) in S. For the Lopatinskii determinant, by [561 Lemma 5.3] and the 
continuity of the Lopatinskii determinant, we can obtain (17.31) and (|7.2p for the invertible 
points and roots of the Lopatinskii determinant respectively. Therefore the argument in 
Section [6] can be applied, and we can derive the result in |56] . 
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